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-DEPENDENT  AMPLITUDE-DISTANCE  CURVE  FOR  P-WAVES  FROM  87°  TO  110° 


JAMES  H.  ANSELL 

Summary 

This  work  is  an  attempt  to  clarify  the  nature  of  the  amplitude- 
distance  curve  for  P-vaves  between  87°  and  110°,  using  thr>  spectral 
amplitudes  of  earthquakes  in  the  Indonesian  region  recorded  at  the  Swedish 
and  Finnish  seismograph  stations.  At  the  present  time  the  results  are 
inconclusive,  because  even  after  allowing  for  station  and  source  terms 
there  is  a  large  unexpected  scatter. 

Introduction 

The  amplitude  decrease  of  P-waves  which  are  incident  at  the  surface 
beyond  95°  and  which  have  passed  through  the  deeper  part  of  the  mantle  has 
been  observed  for  a  long  time.  As  the  quality  and  distribution  of 
seismographs  have  increased,  it  has  been  possible  to  make  more  detailed 
studies.  Gutenberg  (i960)  with  limited  data  demonstrated  the  frequency 
dependence  of  the  amplitude  decrease,  while  Sacks  (1966)  with  better  data 
illustrated  the  effect  very  clearly.  These  results,  however,  are  more 
qualitative  than  quantitative  and  cannot  be  used  to  test  hypotheses  on  the 
nature  of  the  core^iantle  boundary  region  which  produces  the  shadowing 
effect.  Also  these  studies  use  a  few  earthquakes  and  do  not  consider  the 
effects  of  station  geology  on  the  results.  More  recently,  Alexander  and 
Phinney  (1966)  have  worked  with  long-period  waves  in  the  shaacrw  region, 
but  their  data  has  large  scatter,  they  do  not  consider  station  effects  and 
they  do  not  combine  data  from  different  earthquakes. 

Recently,  Carpenter,  Marshall  and  Douglas  (l£wi/  and  Cleary  (1967) 
have  worked  on  the  amplitude-distance  curve  between  30°  and  102°  and  have 
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used  a  joint  analysis  method  described  by  Carpenter  et  al.  This  method 
allovs  the  combination  of  different  earthquakes,  finds  corrections  for  the 
station  effect  and  produces  an  amplitude-distance  curve  independent  of 
earthquakes  and  stations.  These  amplitude-distance  curves  are  valid  for 
short-period  vertical-component  records  of  about  1  aec  period,  but  there 
is  some  disagreement  between  them  at  distances  beyond  90°  probably  because 
of  the  different  methods  of  measuring  amplitudes.  Both  authors  have  little 
data  over  95°  and  neither  investigate  the  effect  of  frequency  on  the 
amplitude-distance  curve.  The  present  work  has  been  done  to  attempt  to 
clarify  the  frequency  dependence  of  the  amplitude-distance  curve  using  the 
joint  analysis  technique  to  combine  data  from  many  earthquakes  and  many 
stations. 

The  observational  material  used  in  the  present  investigation 
consists  of  short-period  vertical-component  records  of  P-waves  from  the 
network  of  Swedish  and  Finnish  stations  for  a  number  of  earthquakes  in  the 
Indonesian  archipelago. 

Analytical  method 

At  teleaeiamic  distances  we  can  express  the  frequency-dependent 
amplitude  A  of  the  body  waves  in  the  form 

A  ■  B.R.S  (1) 

where  B  is  the  frequency-dependent  source  function  which  includes  the 
effect  of  the  crust  and  upper  mantle  at  the  source,  R  is  the  transmission 
coefficient  for  passage  through  the  mantle  which  includes  the  effects  of 
reflection  and  diffraction  by  the  core-mantle  boundary  region  (if  the  wave 
concerned  is  affected  by  this  region),  the  effect  of  transmission  at  any 
boundaries  and  the  effect  of  geometrical  spreading  of  the  waves,  and 
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finally  S  is  the  receive!  function  which  includes  the  effects  of  the 
station  seismographic  response  curve  and  the  crust  and  the  upper  mantle 
below  the  station.  Each  of  B,  R  and  S  also  includes  the  effect  of  the 
anelastic  dissipation  and  scattering  by  inhomogeneities  in  the  regions 
concerned.  B  and  S  vary  with  azimuth  euad  also  with  angle  of  incidence  on 
the  surface.  To  the  extent  that  the  lower  mantle  is  inhomogeneous  R  is 
dependent  on  the  particular  path  through  the  mantle. 

In  our  problem  then,  we  have  selected  the  stations  and  earthquakes 
such  that  we  make  the  assumption  that  B  and  S  vary  little  over  the  small 
variation  of  azimuthal  angles  and  small  variation  of  angles  of  incidence 
involved  (this  assumption  may  not  be  valid!).  R  will  apply  to  the  mantle 
between  Indonesia  and  Fennoscandia  and  to  the  core-mantle  boundary  region 
under  Central  Asia.  All  of  A,  B,  R  and  S  are  frequency-dependent  and 
complex. 

If  we  use  the  base  ten  logarithms  of  these  quantities  then  we  have 

a  ■  b  +  r  +  8  (2a) 

where  a  «  log10|A|t  b  ■  log10)B|  etc.  (|a|  is  the  amplitude  of  the  complex 
A)  and 

phase  (A)  «  phase  (B)  +  phase  (R)  +  phase  (S)  (2b) 

For  any  particular  measurement  of  a  w?  have 

e»b+r+s+e  (3) 

where  e  is  an  error  term  which  includes  the  inaccuracies  of  measurement  of 
a  and  the  effect  of  inadequacies  of  the  model  we  have  set  up.  This 
formulation  is  the  same  as  derived  by  Carpenter  et  al.  (1967)  except  that 
our  £  is  the  log  of  the  spectral  amplitude  and  not  log^g  (^)  and  in  our 
model  the  azimuths  and  angles  of  incidence  are  very  limited  in  range. 
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Following  Carpenter  et  al.  we  find  that  if  we  make  a  number  of 

observations  of  ^  at  a  number  of  stations  for  a  number  of  earthquakes  we 

can  obtain  estimates  of  b,  r  and  s.  If  we  denote  by  subscript  i  the 

particular  earthquake  considered  then  b^  is  the  source  term  of  the  i 

earthquake.  Similarly,  if  we  denote  by  subscript  j  the  particular  station 

considered  then  s.  is  the  station  (crustal  +  seismograph)  function  for  the 
J 

Jthstation.  Finally,  if  we  divide  the  distance  range  into  intervals  over 
which  the  amplitude-distance  curve  is  assumed  constant  and  we  denote  by  k 
the  kth  such  interval,  then  r.  is  the  mantle  transfer  function  for  this 
distance  range.  If  some  estimate  rg  of  the  amplitude-distance  curve  is 
available,  then  we  can  subtract  r  from  both  sides  of  equation  (3)  and  r. 

6  K 

is  considered  as  the  actual  difference  of  r  and  r  .  When  r  is  a  reasonable 

e  e 

approximation  to  r,  the  constancy  of  over  a  distance  interval  is  a  less 
imposing  condition  and  yet  we  retain  the  flexibility  of  the  histogram 
representation. 

Thus  if  earthquake  i  is  observed  at  station  j  and  the  separation  of 
the  two  is  in  distance  range  k,  the  observed  amplitude  a^k  may  be 
expressed  in  the  form 


a. ..  *  b.  +  r.  +  s.  +  e. .. 
ljk  1  k  j  ljk 


(U) 


For  Nre  observations  of  a^^  from  Nep  epicentres  at  some  or  all  of  Nst 
stations  using  Npa  distance  ranges,  we  have  a  set  of  Nre  linear  equations 
for  a-j^.  We  have  Nep  unknowns  b^ ,  Npa  unknowns  r^,  Nst  unknowns  Sj  and 
Nre  unknowns  e..t.  If  we  remove  from  b.,  r  and  s>  their  respective 
averages  so  that 


C  ■  b.  +  r.  +  s. 
i  k  j 


b. 

rk  -7k 


b.' 

i 


8.  -  S.  »  S.' 


(5) 
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then  the  nev  b^',  r^'  and  s^*  averaged  over  i,  k  and  j  respectively  are 
zero.  Hence  ve  have  Nre  +  3  equations: 

Hre  equations 


aijk  *  C  +  V  *  V  *  “o'  +  cijk 


and  3  equations 


(6) 


Z  b . '  ■  0 ,  E  rJ  «  0  and  I  s . '  =0 
i  1  k  0  J 

We  can  henceforth  drop  the  primes  on  b.,  r  and  a., 

1  K  J 

We  can  represent  these  equations  in  the  matrix  formulation 

P  =  QX  +  E 


(7) 


where  P  is  the  row  vector  of  a^k  in  some  order,  E  is  the  error  row 
vector  of  in  the  same  order  as  a^^,  X  is  the  column  vector  (C,^, 

V**  bNep'  rl*  r2***  rNpa*  Bl*  B2**,BNst^  ^  Q  ia  the  mfttrix  of 
indicator  variables  such  that  if  the  nth  element  of  P  is  a. then  the 

ljk 

nth  row  of  Q  multiplied  by  X  gives  C  +  b.  +  r.  +  s.  and  the  last  three 

1  K  J 

rows  of  Q  when  multiplied  by  X  give  equations  (6). 

It  is  possible  to  solve  this  matrix  equation  by  the  least  squares 
method  to  minimise  |e|  and  get  an  estimate  of  X  and  hence  of  C  and  of  K, 

rk  and  v 


The  least  squares  estimate  for  X  is  given  hr 

x  -  (qtq)  1(QTP) 

-1 


(8) 


T  ip  T 

where  Q  is  the  transpose  of  Q  and  (Q  Q)  is  the  inverse  of  Q  Q,  i.e. 

T  t 

(Q  Q)  (Q  Q)  ■  I,  the  identity  matrix.  Problems  may  arise  with 
•  T  T 

calculation  of  (Q  Q)  (Q  p)  and  these  problems  are  discussed  by 


Anderssen  (1969),  In  the  present  work  straightforward  matrix  inversion 
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T  T  T 

was  used  to  fora  (Q  Q)  and  the  difference  (Q  Q)  (qaq)  -I  was  used  as  a 

T 

guide  to  the  accuracy  of  the  inversion  of  Q  Q.  Since  Q  is  composed  of 

T 

integer  indicator  variables,  Q  Q  can  be  calculated  exactly  and  is  not 
affected  by  computational  rounding  errors.  If  there  are  a  sufficient 
number  Ni  of  linearly  independent  equations  (M,  then  a  solution  X  of 
equation  (7)  can  be  found.  Mi  +  3  must  be  greater  than  1  +  Nep  +  Npa  +  Nst, 
and  the  greater  Ni  the  better  the  statistical  estimate  of  X. 

Observational  material 

The  stations  used  are  those  of  the  high  quality  Svedish  and  Finnish 
networks  situated  on  the  relatively  homogeneous  Fennoscandian  shield 
(table  l).  The  earthquakes  used  occurred  in  the  Indonesian  area  between 
the  beginning  of  1963  and  the  end  of  1968  (table  2).  See  also  figure  1. 

The  particular  earthquakes  selected  were  such  that  the  signal-to-noise 
ratio  was  generally  good,  the  amplitude  of  the  signal  was  sufficient  to 
make  further  analysis  worthwhile  and  the  energy  of  the  signal  was 
concentrated  near  the  onset.  Any  selection  of  the  data  will  affect  the 
final  result  as  the  criteria  used  are  subjective.  If,  for  example,  a  record 
is  rejected  because  of  low  signal-to-noise  ratio  -  then  it  may  be  that  the 
noise  level  is  high  or  that  the  'amplitude  level  is  low.  However,  some 
selection  must  be  made  and  the  criteria  used  seem  reasonable. 

The  stations  and  earthquakes  are  related  so  that  for  any  one  earth¬ 
quake  the  stations  cover  an  azimuthal  range  of  less  than  10°  and  that  for 
one  station  the  earthquakes  cover  a  back  azimuthal  range  of  less  than 
about  20°  (cf  figure  2).  For  the  range  90°-110°  epicentfal  distance,  the 
angle  of  approach  of  the  seismic  P-wave  changes  little.  So  for  each  earth¬ 
quake  the  station  net  covers  a  small  solid  angle  and  for  each  station  the 
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earthquake  epicentre*  cover  a  mall  solid  angle.  A*  ve  shall  see  later, 
these  conditions  should  make  the  joint  analysis  method  suitablo  for 
analysing  the  data. 

The  eleven  Swedish  and  Finnish  stations  originally  chosen  are  given 
in  table  1.  Of  these  SOD  was  later  rejected,  because  of  the  nonstability 
of  its  amplification  curve,  and  UDD,  which  because  of  its  later 
construction  recorded  only  four  of  the  earthquakes  (two  on  the  earlier 
Orcnet  instrusent  and  two  on  the  later  installed  Benioff). 

Sixteen  earthquakes  were  initially  selected,  listed  in  table  2. 
Eleven  of  these  earthquakes  lie  in  a  narrow  back  azimuthal  range  from 
Scandinavia  and  the  other  five  are  outside  this  bond.  The  latter  five  are 
treated  as  suspect,  as  the  station  terms  may  vary  too  much  with  large 
changes  in  back  azimuth.  Of  the  original  1?6  possible  records,  109  were 
selected  and  digitised.  Ik  records  were  not  available,  20  were  at  too 
short  epi central  distances  and  33  were  rejected  because  the  signal-to- 
noise  ratio  was  too  low  or  the  record  amplitude  was  not  large  enough  to 
make  Fourier  analysis  worthwhile*  Figure  3  shows  a  typical  record. 

For  each  earthquake  the  epicentral  distance,  azianith  and  back 
azimuth  to  the  stations  of  the  net  were  calculated.  The  epicentral 
distances  were  corrected  for  depth  of  focus  using  the  results  of 
Buchbinder  (1968).  These  corrections  are  such  that  all  the  earthquakes 
can  be  considered  as  surface  focus  events  with  regard  to  the  amplitude- 
distance  curve. 

Experimental  method 

For  the  records  from  each  earthquake  a  suitable  record  length  was 
chosen,  either  20,  30  or  kO  sec,  and  this  length  was  such  that  the  main 
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part  of  the  energy  was  in  the  earlier  portion  of  the  record  and  at  the 

end  of  the  record  the  amplitude  was  much  smaller  or  reduced  to  near  noise 

level*  This  selection  of  record  lengths  should  minimise  the  effects  of 

truncation  of  the  record*  The  start  of  the  record  was  taken  just  bef  -e 

the  apparent  onset  of  the  arriving  P-vave. 

The  records  were  photographically  enlarged  four  or  five  times*  Then 

the  top  and  bottom  of  the  trace  were  digitised  on  a  DMac  pen  follower  and 

the  data  were  converted  to  cards.  They  were  then  interpolated  to  the 

desired  interpolation  interval:  sec  for  20  sec  records,  Sec  for 

UO 

30  sec  records  and  sec  for  bO  sec  records.  The  average  of  the  two 
traces  was  taken  and  the  Fourier  transform  of  the  average  computed  in  the 
form  of  amplitude  and  phase  spectra.  The  theory  of  the  spectral  analysis 
of  digitised  seismic  data  is  well  covered  by  Huang  (1966).  If  the  seismic 
trace  is  the  time  function  f(t),  then  the  computed  Fourier  spectrum  is 
given  by  F(vq)  where 


F(vn) 


1 

m 


m-1 

I  f(kdt)  e 

k-0 


2irinkAt 

““"t 


(9) 


T  is  the  length  of  the  record,  vQ  is  the  nth  frequency  in  cycles  per  sec 
and  vn  ■  ^  where  n  runs  from  0  to  m,  m  is  the  number  of  digitised  points, 
At  is  the  digitising  interval  and  mdt  ■  T,  and  finally  for  most  efficient 
conputation  m  is  a  power  of  2  (in  our  case  m  ■  256  or  512).  We  avoid 
aliasing  by  using  a  digitising  interval  sufficiently  small  so  that  the 
amplitude  spectra  are  neglible  for  frequencies  above  the  folding 
frequency  ~  . 

Various  methods  of  windowing  were  considered  but  none  was  applied 
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as  none  seemed  suitable.  Using  longer  record  lengths  vith  low  cut-off 
amplitudes  should  minimise  the  effect  of  truncating  the  records. 

One  record  of  average  quality  was  photographically  enlarged  and 
digitised  separately  three  times.  The  agreement  between  the  three 
amplitude  spectra  is  very  good  as  shown  in  figure  4.  The  maximum 
variation  throughout  most  of  the  frequency  range  was  0.25  units  compared 
with  the  maximum  amplitude  of  4  units.  For  the  larger  amplitude 
components  the  difference  is  less  than  8  %,  For  better  quality  records 
the  agreement  should  be  better  and  the  opposite  for  poorer  quality  records. 

The  amplitude  spectra  have  been  smoothed  using  a  3-point  smoothing 
with  ■£,  weighting  for  the  20  sec  records,  a  3-point  smoothing  with 

—  weighting  for  the  30  sec  records  and  a  5-point  smoothing  with 

it*  V*  V  I'  weighting  for  40  sec  records.  This  smoothing  should  improve 
the  consistency  of  the  results  and  also  make  the  comparison  of  the  spectra 
frem  different  length  records  more  meaningful.  This  smoothing  is  not 
windowing  but  an  attempt  to  smooth  the  insignificant  fluctuations  in  the 
amplitude  spectra. 

All  the  spectra  obtained  are  divided  by  the  instrument  magnification 
factor  at  1  sec,  and  so  the  amplification  curves  are  normalised  at  1  sec. 
Also  the  inverse  of  the  magnification  factors  for  the  photo  enlargement 
is  applied  such  that  the  amplitude  is  in  units  of  0.1  microns  and  after 
we  have  taken  the  log^  of  the  amplitude  spectra,  we  add  one  to  the 
results,  i.e.  the  log  (amplitude)  of  the  spectra  is  such  that  the 
amplitude  is  measured  in  0.01  microns.  PcP  is  always  included  in  the 
pulse,  and  if  the  earthquake  is  shallow,  pP  is  included  in  the  record.  If 
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the  earthquake  is  deeper,  pP  either  does  not  affect  the  record  or  is 
small  and  appears  near  the  end  of  the  record. 

Computations  and  results 

In  the  experimental  work  we  find  estimates  of  a^  for  various 

earthquakes  and  stations.  Then  we  apply  the  analytical  method  of  joint 

analysis  to  estimate  the  amplitude-distance  curve  and  the  station  terms. 

A  computer  program  to  solve  equation  (7)  and  find  X  in  the  form  of 

equation  (8)  has  been  developed.  The  program  calculates  the  station  and 

source  terms  s.  and  b.  and  also  the  amplitude-distance  curve  using  a 
J  1 

histogram  of  2°  intervals.  See  Appendix. 

The  results  are  very  poor.  In  figure  5  we  show  a  plot  of  the  raw 
data  from  which  is  subtracted  the  appropriate  source  and  station  terms  and 
the  constant  introduced  in  equation  (5).  Even  though  the  station  and 
source  terms  are  allowed  for,  the  scatter  is  very  high  and  certainly  no 
amplitude  decrease  is  seen  beyond  90°  -  as  would  be  expected.  This 
behaviour  seems  to  come  from  the  data  and  not  from  the  iuversion  program. 
So  far  no  explanation  has  been  found  for  the  anomalous  behaviour  of  the 
results.  The  data  has  been  divided  into  smaller  groups  of  earthquakes 
with  narrower  azimuthal  and  distance  ranges  but  there  is  no  substantial 
improvement  in  the  results.  It  is  possible  that  the  model  we  have  set  up 
is  based  on  invalid  assumptions  on  the  nature  of  the  source  and  station 
functions. 

As  an  example  we  present  the  amplitude  data  for  1  sec  period  for 


all  stations  which  recorded  the  earthquake  on  the  29  July,  1968.  We  list 
the  stations  in  order  of  azimuth  (epicenter  to  station)  and  epicentre! 
distance.  It  is  obvious  that  no  clear  pattern  emerges. 
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Station 

Azimuth 

log10(Amp)  at  1  sec 

KLS 

329?1 

-  0.1664 

GOT 

331.5 

-  0.1292 

NUR 

331.7 

-  0.1053 

UPP 

332.2 

+  0.36l 4 

UDD 

334.0 

-  0.4820 

KJN 

334.8 

-  0.0336 

UME 

335.5 

-  0.3136 

SKA 

336.8 

-  0.4511 

KIR 

339.3 

+  0.2392 

Station 

Distance 
(reduced 
to  zero 
focus ) 

lOe10(Amp)  at  1  sec 

KJN 

97°1 

-  0.0336 

KIR 

98.7 

+  0.2392 

NUR 

99.4 

-  0.1053 

UME 

100.3 

-  0.3136 

UPP 

102.9 

♦  0.361 4 

SKA 

103.6 

-  0.4511 

UDD 

104.3 

-  0.4820 

KLS 

105.3 

-  0.1664 

GOT 

106.5 

-  0.1292 

Again  for  the  earthquake  on  the  1$  July,  196!/,  we  have  the  following 

results  at  1  cycle  per  sec  frequency* 

(Obviously  the  data  is  not  accurate 

to  four  decimal  places! ) 

Station 

Azimuth 

log1Q(Amp)  at  1  sec 

KLS 

327?8 

0.1578 

GOT 

330.0 

0.2462 

UPP 

331.0 

0.2488 

UME 

334.4 

0.1685 
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SKA 

335.1* 

-  O.U373 

KIR 

338.3 

0.7596 

Station 

Distance 

log10(Amp)  at 

KIR 

90?1 

0.7596 

UME 

91.3 

0.1665 

UPf 

93.6 

0.21*88 

SKA 

91*. 7 

-  0.1*373 

KLS 

95.7 

0.1578 

GOT 

97.0 

0.21*62 

In  the  first  example  the  back  azimuths  vary  from  61°  to  75°  and  in  the 
second  example  from  67°  to  7^°.  The  change  in  back  azimuths  between  the 
two  examples  is  about  5*5°  for  each  station. 

The  Fourier  spectral  program  from  seismogram  to  amplitude  spectrum 
has  been  checked  against  an  independent  program.  The  spectral  estimates 
seem  therefore  to  be  valid. 

The  problem  remains  -  which  of  the  assumptions  ve  lave  made  is  not 
valid?  Possible  it  is  that  the  source  function  can  vary  very  rapidly  over 
very  small  azimuthal  angles.  In  both  the  above  examples  the  smallest  and 
the  largest  amplitudes  are  next  to  each  other  in  the  distribution  of 
azimuth.  (This  effect  has  not  been  checked  on  other  data  sets).  If  the 
source  spectrum  does  vary  so  much  with  such  small  angles,  then  spectral 
analysis  of  short-period  P-vaves  from  earthquakes  could  only  be  done  on 
a  statistical  basis.  Explosions  provide  much  more  symmetrical  sources, 
but  their  limited  distribution  prohibits  their  application  to  our  present 
problem* 
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Table  1 


Btationa  uaed  (aee  also  figure  1) 


Station  Code 

Station  Name 

Location 

Latitude 

Longitude 

UPP 

Uppi&la 

59°86n 

17°63E 

UME 

Utaei 

63.82 

20.24 

KLS 

Karlskrona 

56.17 

15.59 

GOT 

Odtebore 

57.70 

11.98 

UDD 

Uddeholm 

60.90 

13.61 

SKA 

Skalatugan 

63.58 

12.28 

KIR 

Kiruna 

67.84 

20.42 

KEV 

Kevo 

69.76 

27.01 

SOD 

SodankylA 

67.37 

26.63 

NUR 

Nuntij&rvi 

60.51 

24.65 

KJH 

Kajaani 

64.10 

27.70 
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Table  2 


Earthquakes  used  (see  also  figure  1) 


Data  from  USCGS 


Date 

Origin  time 
(GMT) 

h  m  s 

Epicentre 

Latitude  Longitude 

Depth 

(km) 

Magnitude 

(m) 

(UPP.KIR) 

26.2.1963 

20.14.08.7 

-7?5 

l46?2E 

171 

7.7 

7.4.1963 

22.36.03.4 

-4.9 

103.2 

72 

6.7 

21. 3. 196 4 

03.42.19.6 

-6.4 

127.9 

367 

6.6 

28.3.1964 

11.30.09.8 

0.5 

122.3 

140 

6.2 

8.7.1964 

11.55.39.1 

-5.5 

129.8 

165 

7.1 

18.10.1964 

12.32.24.1 

-7.0 

124.0 

574 

7.0 

29.4.1965 

15.48.57.1 

-5.6 

110.2 

504 

6.3 

15.7.1965 

18.33.29.9 

7.7 

123.8 

588 

6.5 

20.7.1965 

13.18.27.4 

7.5 

124.3 

45 

5.9 

20.8.1965 

05.54.50.0 

-5.7 

128.6 

326 

6.7 

21.8.1966 

05.00.26.6 

8.5 

126.7 

67 

6.2 

21.5.1967 

18.45.11.7 

-1.0 

101.5 

173 

7.0 

26.8.1967 

00.36.42.1 

12.2 

140.7 

33 

6.6 

24.5.1968 

15.43.54.2 

-6.8 

116.9 

605 

6.3 

29.7.1968 

23.52.15.0 

-0.2 

133.4 

12 

6.7 

27.9.1968 

03.58.55.1 

-6.8 

129.1 

127 

6.9 

Figure  captions 
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Fig*  1*  Mercator  projection  of  area  of  interest  (shoving  Fennoscandian 
stations  and  Indonesian  epicentres  used). 

Fig*  2.  Cross-section  of  the  earth  shoving  diagrammatically  the  ray  paths 
from  Indonesia  to  Fennoscandia  (the  diagram  is  merely  suggestive 
and  not  to  scale)* 

Fig*  3*  Short -period  vertical-component  P-vave  recorded  at  Ume&  from  the 
Banda  Sea  earthquake  of  21  March,  196b. 

Fig.  b.  The  effect  on  the  spectral  amplitude  of  treating  the  record  at 
Kiruna  from  the  earthquake  of  21  March,  196b,  three  times  as  a 
separate  unit* 

Fig.  9*  Plot  of  all  rav  data  for  1  cps  vith  station  and  epicentre  terms 
removed  (shovs  failure  of  method  to  give  the  expected  result  and 
also  shovs  large  scatter). 


stations 


FIG.  1 


FIG.  3 


vvwwo^n. 


■■  i. 


LOG..  AMPLITUDE  0.01 


fr#q.  1  c.p.s. 


i? 

CO 


DISTANCE 


Appendix 


Program  SOLVE  is  used  to  invert  the  system  of  linear  equations 
described  in  the  section  "Analytical  method".  The  program  is  not 
particularly  efficient  but  is  effective.  Subroutine  DIST  1  puts  the 
earthquake-station  pairs  into  their  distance  ranges.  Subroutine  SHUFFLE 
is  used  to  vary  the  input  without  changing  the  coding  of  the  original 
seismograms  -  earthquake  3  at  station  6  has  code  1316.  Subroutine  SELECT 
selects  the  frequencies  required  from  the  20  frequencies  of  the  input. 
Subroutine  DIST  2  is  here  only  formal  but  may  be  use^  to  remove  any  prior 
estimate  of  station,  distance  or  source  terms. 


r>  o  n  .  r<  n  n 


'''JOfl,  1C4UC  «A\b*-LL  .  3 
*DCVAND, 2SOOO 

^TN^.D.p.L  »X 

PPOr-DAM  S^LV/p 
r  ALL  WP9< 


NOT  REPRODUCIBLE 


Fa  se  1A 


"ALL  EXIT 

pnd 

e;i IPPmjT  J  Nc  voo< 

THIS  ORCr-ffAV  rc.jLO  BP  MORE  EFFICIENT 

MAXIMUM  D I VPNS I  'JiN  b  ll'J  KfcCOKOS  ♦  15  FREO’JEC  ItS  AND 

NST  STmTiONS  «  NkA  TLr\b  AND  NL^  cc  i  CLw I K«_b  WHtKt 

NST+nPA+ulP  iS  LEbb  THmN  OR  tv*U«u  0  49 
DIMENSION  A(  113,50)  ,5<  113,15)  .DISTANCE)  1 10) «FKQ<20) «ANUK(60) 

D  |  ^ENSIGN  XI  50.1 5)  .AS  I  50. 50.  AT)  125,20) 

OT VPNS  I  ON  IDR( 113) «ATCOL( 11?) .ITITLE(6) 

POUT.VALLNCi-  (  AS.  AT  i 

C^VMfN/lC/AA( So,  50  ) 

RFAD  991  ,NST ,NrP«NPA.NPP«NPR 
'JQl  porvat(5I5) 

ORINT  992  .NST.NEP.NPA.NRE.NFP 

9*2  FORMAT) lx, 15,*  STATIONS.  *,I5,*  EPICENTRES,  *  *  I  5 . *  PARAMETERS, 

1  *,!5,*  RcCORPS«  ♦.15,*  PREOUFNCIES.  *) 

NC0L*NST+NPP+NDA+1«  NK0W*NRt+3 
READ  1 f  (PPO( I »,i«l.NFR) 

1  FOPM*r< 15F5.2) 

KST»1+NSTSKPA3 1+NST+NPAS JST  a2S JPA« 2+NSTi JEP*2+NST+NPA 

<cp*^roL 

I TITLC< 1 j.flMFPFO. 

DO  505,1*3,6 

3r»5  I T I  Tl  :  ( I ) “PH 

OUT  A  (  I  «  J )  «0.0 
On  2,r»l,NR0W 
Do  3, j*l ,NCOL 
A ( T,J)“O.D 

3  CONTINUE 

2  CONTTMUF 

PUT  IN  VALUFS  Oc  FIRST  COLUMN  0  AU,J)  AND  BOTTOM  THREE  ROWS  OF 
A  AND  B 
Do  4,J«1.N0P 

4  A(  I  ,H«I  ,0 

NRct»N9P*1 

D0~5, I-JST.KST 

5  A(NRPT,n*1.0 

NRFT.NPFT+? 

Do  6,1 “JEP.NCOL 

6  A ( NRCT » I ) ■ 1 ,0 
nrft«nrft-i 

DO  9001 , I » JDA.XPA 
oooi  A(N9ET,H«1.0 
NRPT«N9F+1 
DO  7,J«1,NFR 


DP  «,  I*NRETtNPPV 
B(  I  •  J)*n*o 
8  CONTINUE 

7  CONTINUE 

ANUMI I »  GIVES  THE  NUMBER  OF  NON-ZERO  ELEMENTS 
POINT  12 

12  PPRMAT(4X.*INPUT  DATA*) 

PRINT  14 

14  FORMAT ( *  *1 

PRINT  15*{PROf  D.I-l.NFR) 

15  PORMAT(Rx.*FRPOUFNCY  *  * 15F7#  2 ) 

POINT  14 

14  FPRM^nlX.*  CPI  ST  XT  DISTANCE*) 

D«  0,1  «2*NCr'L 
0  ANIJM  (  1  1  «o.r» 

ANUM{ 1 )«PLOATP(NRP) 


C  rpad  in  DATA 

OP  100* I *1 *NRp 
READ  10.MEPIC.MSTAT*DISTA 
10  PPRMAT (1X*2I2*F6*1> 

DI STANCE  ( 1 ) *D I STA 
MoPir«MFPIC-10  5MSTAT«MSTAT-1^ 

Mr>UMP«MSTAT 

CALL  CHUFFLP(MCPK,NPP*MSTAT*NST) 

M«M<ll\T*l  5N*MF°  I  C*<PA 
A  « | ,M)«1 .0 
A(I *M)«1 .« 

ANUM(  j|)«ANUM{M)4.1,0 
ANUM(N)«ANU^(N)+1.0 
CALL  OISTKDISTA.NPA.KDIST) 

KA»K0TST+l4-NSl 
IPPII )«<A 
A ( I »KA ) ■ 1 #0 
ANUM(XA)»ANIJ^MCAU1,C. 

OPAO  71*  (  FR(J{  <)  ,  X»1  *2C) 

28  FnRVAT(5F12.4) 

CALL  SFLPCT ( PPO»NeR ) 

C  PRO  IS  HERP  the  input  amplitude  DATA 

PRINT  24*MPPIC.MSTAT.DISTAt<FRQ(<)*K«l,NFR> 

74  PnRMAT IlX»2I5*Ffl*l»llFlO«61 
Z-DISTANCPII ) 

C  SUBTRACT  FITTEO  CjRvE*  SPECIFIED  IN  DIST2,  AND 
CALL  0IST2(FR0,Z,NFR,MSTAT,NST1 
DP  0081 .<*1 *NPR 
oo«l  BC I ,K )«FRO(<) 
inn  CPNTTNIJF 

C  HAVE  MOW  PED  DATm  INTO  B  AND  PARAMETERS  INTO  A 
C  FITTED  CURVE  ANP  INST  EFPFCTS  FRQMB 


Page  2* 

N  EACH  COLUMN  OF  A 


INST  EFFECT 

AND  SUBTRACTED 


op  f  nt  loni 

'0"1  c^PMAT  (*  MATRIX  A*1 

no  ions, |=itN°nw 

DPINT  lno2*( A( I  ,J1 *J=l*NCOLi  P85e  3<fc’ 

irn^  rr>»jTj^i|C 

’no2  Co°MAT  {  lX»^^c?»r' * 

0  0  I N  T  u 
PRINT  1004 

’ 004  FORMAT (  *  MATRIX  B*1 
no  inns  , t  =  !*npow 
OP  I  NT  1 nOS  * { B( I  «J)*J=l*NFPi 
inn*  CnNTI  N'JF 
lnnfc  FORMAT  {  lX*lSr0«3) 
op  TNT  14 

C  WF  NOW  FCPM  ATA  AND  ATB  SATA  1£  1 NAA  AND  ATB  IS  IN  X 

D"  1(  ♦  .  IM*NCr'L 
no  to*,j  =  i*t 
A  A  (  T  *  n«o.n 

no  USt<*lfW^W 

1  f»  AA(  !  ,J1=AA(I  ♦Jl+AtKf !  >»A(<*J) 

AS ( T ♦ J I  *  A A { T »  J> 
in*  OoMTfN(jF 
I.  *  I  - 1 

On  in^,jri,L 
AAfJ, TlsAAC! ,J) 

’07  AC( j,T!*AA( J,|  I 
104  oontT'MUc 

no  110  *1=1 *N^OL 

noi i i ,J=1,MCP 
X( I ,Jl*O.M 
00  11?*K*1*  MRO'V 

112  X(!»Jl*X(!tJl+A(K»n*B(K»J) 

111  r^NTJNiJF 
1  1 A  fOMJT  V'JF 
ootmt  ln]o 

iftlo  FOPVATC*  MATRIX  ATA* I 


loll 


oqo2 


no  inil,t*l,NCOl 

OP INT  100?* < AA( 1 ,J) *J»1*NC0L> 

OOAfTINUF 

DFTPRM*n,o 

NMAX  =  ^(? 

CALL  MAT!NV(NCOL*X*NFR*OETERM,NNAX ) 

DPINT  ?PV2*DETPPV 

format ( IX *«OFTFPM  =»*E12*4) 

r.pF-o.o 

no  ?r>oi*i  =  l*Nroi 
on  ?rr)?, J»T  *NCA». 

VARso.<) 


NOT  REPRODUCIBLE 


o  o 


Do  2f  d**k* ! ,Nrr>L 

’  VAR*''*R+AA(  f  ,<)*AS( K* J) 
IFCI.iQ.J)  VAR«VAR-1.0 
VA9«VAR#VAR 

2002  IFCVAR.GT.GRF)  GRE=VAR 
2001  CONTI  NUE 
ORT^T  14 
r,oc,SORT':|r-DC) 


NOT  REPRODUCIBLE 


Tage  4**- 


ooiMT  ?oo*,r.Rr 

20  5  FORMA T ( *  GRFAT-ST  FRROR  IN  PRODUCT  OF  MATRIX  AND  ITS  INVERSE  * 


*. 


1E12.4/) 

DO  120*1*1 »NROW 
DO  121  ,J«1,NFP 
O'IM.0,0 

DO  1??,  K= 1 *NrOL 
122  DUM  »DUM+A (I«K)*X(K*J) 

AT( T  ,J)*DUM-B<  I.J) 

B( I  * J ) *AT ( I  * J1»AT ( I ,J> 

ATC I ,J)»XC IORI I  1 » Jl-ATC I , J) 

121  CONTI MUC 
CONTI  M'jF 
DO  A^ft* J*1 *NffP 
DOAo2*!*l,NO* 

*02  ATCot  ( I )  *AT ( I  •  J 1 
ITITI  ~C2)*AH 
CALL  :NC0DF( ITITLF) 

CALL  FMTS(P) 

CALL  FMTKJ.l) 

CALL  GRAPH1( DISTANCE, ATCOL .-NRE.3H7X8 *4HAuT0*I TITLE *iOHDI STANCE i - , 
15HAMP.#,526060fco6O606O6oB} 

*00  CONTI 4UE 

X  CONTAINS  THC  SOLUTIONS  AND  9  CONTAINS  THE  ERROR  SQUARED 
WF  FIND  NOW  THE  STANDARD  DEVIATIONS 
SOANl *SORTF( ANUM( 1 ) -1 . 0 ) 

SDAN2-S0RTFC ANUM1  1 >-NPA*l. 0-1.0) 

ANNN*  ANUM ( 1 1-1.0* | NPA+NSl+NEP+l I 
ANNN  *MAX1F ( ANNN* 1*0) 

S0AN3«S0PTF( ANNN1 
DO  1  SO  ,  I *1 ,NFR 
AACI, 11*0.0 
Do  151  .  J*1 ,NRR 
1*1  A  A  Cl  •  1)*AA(  I  *  1 )  -f  B  (  J  •  I  ) 

DUM*SORTF(AACI,l> 1 
A  A  (I,!  1*D‘JM/S0AN1 
A A(  I  ,  2)»P'JM/SQAN2 
A A  C I , 3 1 *DUm/SDAN3 
1*0  CONTINUE 

C  EXCEPT  FOR  FIRST  COL  ANUMC I )  IS  NO  IN  COL  MINUS  ONE 
D 0  1*0  ,J«1,NC0L 
<*  J+3 
L*  J*1 

DO  160  *  I  * 1 *NFR 


NOT  REPRODUCIBLE 


=0.n 

HO  k  ,  Mr  J  ,  Mf>p 

Ifl  A  A  (  I  *  < ) = AA (  f  ) *A ( M  »L )  A 

i^n  aa(  i  ,<i=sgrtf(aa<  t»K)/anum(L)  i  Pa*c 


c°  CONT  ’  MIJF 
POTV  14 

noiwr 

’  0  "p SUL  t  S* / ) 

OOTA'1  2r>l 

Pol  POPMATI*  CONSTANT*) 

?r»2  RflRr  AT(  lX,lf>Ffl.4> 

On  204  ,J=1,S 

OPTNT  2o2*(X( 1,1  I »AA( I ,J) ,  I  =  1,NFR) 

0o4  roNTT>IUP 

DOT^T  14 
OOJmT  ■>  r>  «S 

cnRMf r ( 1X,*STATT0N5*> 

’''1  poomaT(1X*T?I 
no  20A  » I  =  2 » XST 
j=t-l  f l  =  T+2 
POINT  207tJ 

PRINT  202,  I  X  <  I,K>,AA(K,L),K*1,NFR) 

’06  continue 

DRTNT  14 
PRINT  2o« 

201  POpMAT( 1X,*0ARAVPTFRS*) 

0O  ?0O  f 1 r JO A  ,<DA 
JM-l-NST  «•  L*T  +  2 

or»iwT?o7t  j 

ORINT  202,<X<  T,<),AA(K,U*K*1#NFR) 

JOO  roMTTMUF 

OR  TNT  1 4 
PRINT  211 

’ll  FORMAT! IX, "EPICENTRES*) 
no  212  ,I*JFP,KFP 
J=  T-l -NST-NOA  *L«I  +  2 
RR1MT207, J 

PRINT  202, (X( I,K>,AA{K,L),K«1,NFR) 

212  CONTINUE 
RETURN 
*mD 

subroutine  DIST10ISTA,NPA,KDISTJ 

X  Is  DISTANCE  TO  RIGHT  OF  FIRST  INTERVALS  IS  NUMBER  OF  FIRST  L0N( 
1  INTERVAL,  y  IS  THE  DISTANCE  BETWEEN  THE  TOP  OF  FIRST  AND  BOTTOM 
1  OF  FIRST  LONG  INTERVAL,  W  IS  WIDTH  OF  SHORT  INTERVAL  ,  LONG  Is  5. 
x*lnl .0  SN«4  SY*4,C  SW*2.0 
7*DI 5  TA-X 

IF(Z.LT,0,0)Z«0,125 


IF(7.0T.YI  r, 0  TO  1 
7*?/V 

<dist*intf<7)+i 

PFTtJPN  P,s*  6Jl 

1  Z*(Z-Y)/S.O 

KOI  ST«INTF(7  )+N 

PPTllFN 

PNO 

SUBPOUT  INF  01 ST7( FRQ.Z  #NFR,MSTAT ,NST ) 

DIMENSION  FPO(l) 
on  to  7 

2  ConTINUP 
Rc-TllO^ 

END 

SUBROUTINE  SHUFFLE (MFP I C,NEP.MSTAT.NST  ) 

OIMFf SION  NUTT( 16) ,NURT( 111 

SORTS  EPICPNTRFS  AND  STATIONS  INTO  NUMBERCD  ORDER 

DATA( (NUTT ( I  I ,1*1,16 )»7,C.O. 0,1 »C,C, 0.0 ,2.3, O.S, 0,4) , ( (NURTI I ) , 1*1. 
l»ll)el»2.3»4,S»6.7.8»R»0.1C) 

MPOjr.NUTT  (MFPfrl 
mST*T*NUPT(MSTAT’ 

PFTIJPN 

PND 

SUBROUTINE  SELECT ( FRO. NFR I 


DIMcnsjon  FPom,FLI20) 

:  SELECTS  THE  FF.EOUENCI ES  REQUIRED  FROM  THE  20  READ 

DATA< (NUM(I), I *1,10 1 *1,2,3,4,8,6,7,8,9,10) 

Do  1  1.1,20 

1  eL(T)«FRO(T) 

DO  7  1.1, NFR 

2  FROm«FL(NUM(  I)  » 

PFTMON 

PND 

SCOPE 

'LOAD, 

(  1,  I 2( )*)R 


UNCLASSIFIED 


DOCUMENT  CONTROL  DATA  R  &  D 


<  ta  //  •  ..f  r,  \  ./  i 


>  n  ,it  tn  •  r  •••.  ’iiii  f  *  r  ■ 
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Thia  Work  Is  an  attempt  to  clarify  the  nature  of  tha  amplltudo-dlatanoa  curve  for 
P-waves  betwean  87  and  110  ,  ualng  the  epectral  amplitudes  of  aarthquakes  In.  the 
Indonesian  ration  recorded  at  the  Swedish  and  Finnish  Belemograph  stations »  At  tha 
prssant  time  the  results  are  Inconclusive,  baoauee  even  after  allowing  for  etatlon 
and  aouroe  terme  there  la  a  large  unexpected  ecattsr. 
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